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ABSTRACT 

o 

• We have carried out near-infrared JHK spectroscopy of a gravitationally lensed 

submillimeter galaxy SMM J14011+0252 at z = 2.565, using OHS and CISCO on the 
Subaru telescope. This object consists of two optical components, Jl and J2, which are 
lensed by the cluster Abell 1835. Jl suffers additional strong lensing by a foreground 



(N 



galaxy at z = 0.25 in the cluster. The rest-frame optical Ha, H/3, and [O n]A3727 
lines are detected in both Jl and J2, and [N n]AA6548, 6583 lines are also detected 
in Jl. A diagnosis of emission-line ratios shows that the excitation source of Jl is 
stellar origin, consistent with previous X-ray observations. The continua of Jl and 
J2 show breaks at rest-frame 4000A, indicating relatively young age. Combined with 
optical photometry, we have carried out model spectrum fitting of J2 and find that 
it is a very young (~50 Myr) galaxy of rather small mass (~ 10 8 Mq) which suffers 
some amount of dust extinction. A new gravitational lensing model is constructed 
to assess both magnification factor and contamination from the lensing galaxy of the 
Q-f component Jl, using HST-F702W image. We have found that Jl suffers strong lensing 

with magnification of ~ 30, and its stellar mass is estimated to be < 1O 9 M . These 

*_ i 

results suggest that SMM J14011+0252 is a major merger system at high redshift that 
undergoes intense star formation, but not a formation site of a giant elliptical. Still 
having plenty of gas, it will transform most of the gas into stars and will evolve into a 
galaxy of < 1O 1O M0. Therefore, this system is possibly an ancestor of a less massive 
galaxy such as a mid-sized elliptical or a spiral at the present. 

Subject headings: galaxies: formation — galaxies: starburst — gravitational lensing 
- galaxies: individual (SMM J14011+0250) 
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1. Introduction 

Since the discovery of a large population of Lyman break galaxies (hereafter LBGs; Steidel 
et al. 2003), a question whether we are seeing the majority of high-z stellar population or more 
stars are hidden in dusty star-forming galaxies has been raised. It is expected that massive galaxy 
formation at high-z may cause rapid chemical evolution and dust enrichment, leading them to a 
dusty phase like ultra-luminous infrared galaxies (hereafter ULIRGs) emitting most of their energy 
in the far-IR. Such objects are supposed to suffer from severe dust extinction, be unable to be 
picked up as LBGs defined by a spectral break in the rest-frame UV and only be detected in the 
far-infrared to submillimeter wavelength. 

Deep submillimeter surveys by SCUBA on JCMT have revealed a large population of dusty 
galaxies at high-z, whose far-IR luminosity exceeds 10 12 L satisfying the criterion of ULIRGs (see 
Blain et al. 2002, and references therein). Their optical to near-IR (NIR) color is often red, and 
some of them display no optical/NIR counterpart even using 8- 10m telescopes (Ivison et al. 2002; 
Fox et al. 2002; Frayer et al. 2004). X-ray observations reveal that some of them contain type-2 
AGNs, but contribution to the total luminosity is not clear (Alexander et al. 2003). Submillimeter 
to radio flux ratios indicate that these objects lie at z > 1 (Lilly et al. 1999; Gear et al. 2000; Hall 
et al. 2001; Smail et al. 2002; Ivison et al. 2002; Fox et al. 2002), and spectroscopic follow-ups 
of radio-selected sources reveal their median redshift to be around 2.4 (Chapman et al. 2003). 
These objects are therefore thought to be an important population to probe the star-formation 
history of the universe. Due to their red color and faintness, their spectroscopic observations are 
time-consuming, so there is only a small sample of spectroscopically observed objects in the NIR 
(Ivison et al. 2000; Frayer et al. 2003; Smail et al. 2003a, b; Tecza et al. 2004; Simpson et al. 
2004). However, spectroscopy in the NIR is important because they cover the rest-frame optical 
wavelength, which provides us important information of their physical status, ionizing photons and 
stellar population. We therefore carried out low-resolution NIR spectroscopic observations of SMM 
J14011+0252, one of the brightest submillimeter galaxies at high-z. 

SMM J14011+0252 is discovered behind the cluster Abell 1835 (z = 0.25) as a weakly lensed 
object (Smail et al. 1998; Barger et al. 1999; Frayer et al. 1999). Intrinsic far-IR luminosity 
of 7 x lO 12 I/0 indicates that this is an ULIRG, and inferred star formation rate (hereafter SFR) is 
fewxl0 3 M Q yr" 1 (Ivison et al. 2000, hereafter 100). It consists of two optical/NIR components; 
Jl and J2. They both lie at the same redshift of z = 2.56, and the separation is 2 / .'l. Both are red 
in the rest-frame optical (R — K = 3.8 for Jl and 3.2 for J2) and Jl is also red in the rest-frame UV 
(U — R = 2.1) while J2 is bluer (U — R = 1.1). Both show only weak Lya emission (100). No other 
emission line is detected in the rest-frame UV, but strong Ha+[N ll]AA6548, 6583 lines exist in the 
NIR spectra (100; Tecza et al. 2004). HST-F702W image reveals a more complex structure 
of Jl; it consists of several bright knots surrounded by an extended diffuse component with an 
ERO component (Jin) elongated toward north (Ivison et al. 2001). Integral field spectroscopy 
in the JHK-bands discovers an extended Ha emission cloud toward Jin (Tecza et al. 2004). A 
continuum break at rest-frame 4000A is found in Jl, which indicates recent (~ 200Myr) massive 
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star-formation activities. 

Radio observations detect CO(3 — 2) (Frayer et al. 1999; Ivison et al. 2001; Downes & 
Solomon 2003) and CO(7 — 6) (Downes & Solomon 2003) lines at z = 2.565, which imply strong 
star-formation activity. The detection of CO(7 — 6) indicates that temperature of the gas cloud 
is 30 — 40K and its density is nn 2 ~ 2000 cm~ 3 (Downes &: Solomon 2003). The CO emission is 
extended over 2'/2x<0 / . / 5 and centered at Jl (Downes & Solomon 2003). A deep X-ray observation 
by Chandra shows no sign of X-ray emission (Fabian et al. 2000), suggesting that this object 
doesn't harbor a luminous AGN. 

In addition to above, a detailed re-analysis of the optical spectrum reveals several absorption 
features at z = 0.25 including Ca H+K (Smail et al, private communication (2003); Downes 
& Solomon 2003). This suggests that the main component of J 1 is a foreground evolved galaxy, 
possibly an elliptical galaxy in Abell 1835, and that the other knots of Jl are strongly-lensed images 
of the submillimeter galaxy at z = 2.565. Downes &: Solomon (2003) estimated the magnification 
factor to be ~ 25 from their CO observation. 

Here, we report results of low-resolution NIR JHK spectroscopy of Jl and J2, using OH- 
airglow suppressor (OHS; Iwamuro et al. 2001) and Cooled Infrared Camera and Spectrograph 
for OHS (CISCO; Motohara et al. 2002) installed at the 8.2m Subaru telescope (lye et al. 2004). 
Thanks to the high sensitivity of OH-airglow suppression technique combined with the 8m aperture, 
high-S/N continuum spectra from rest-frame 3000 to 6500A with several strong emission-lines are 
obtained. 

We describe our observations and data reductions in section 2, display results in section 3, 
construct a new gravitational lensing model and discuss implication to the nature of this object in 
section 4, and summarize them in section 5. Hq = 71 km s _1 Mpc _1 , f^M = 0.27, and Qa = 0.73 
are assumed throughout this paper. 

2. Observations and Data Reductions 

2.1. if-Band Imaging and Spectroscopy 

.ST'-band imaging and if-band spectroscopy were carried out on 2002 February 28, with CISCO 
mounted at the IR-Nasmyth focus of the Subaru telescope. The optical-Nasmyth secondary mirror 
was used, and the pixel scale was O'.'lll pix^ 1 . 

If '-band images were taken with single exposure time of 20s and the telescope was nodded 
after every 3 frames by 10" in a 4-point dither pattern. A total of 12 frames was acquired, for a 
total exposure time of 240s. The seeing size was 0'/45. Before the if'-band imaging, a faint standard 
star FS 27 (Hawarden et al. 2001) was observed for flux calibration. 

The If -band spectra were acquired with / . / 65 slit, providing wavelength resolution of 400 at 
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2.15//m. The position angle of the slit was set at —69° to introduce both Jl and J2 on the slit. 
Single exposure time was 200s, and the telescope was nodded 4"— 6" after every 3 frames for sky 
subtraction. A total of 24 frames was taken for a total exposure time of 4800s. Then, AO star SAO 
120250 was observed for correction of the atmospheric and instrumental transmission curve. 

2.2. i^-Band Imaging and Jii-Band Spectroscopy 

ii-band imaging and JH-band spectroscopy were carried out on 2002 March 1, 2, and May 18, 
with OHS and CISCO mounted at the IR-Nasmyth focus of the Subaru telescope. The opt-Nasmyth 
secondary mirror was used, and the pixel scale was O'.'lll pix -1 . 

.ff-band images were taken with OHS in the imaging mode (see Iwamuro et al. 2001) before 
the spectroscopic observations. Single exposure time was 50s and the telescope was nodded after 
every 3 frames by 10" in the direction of the dark lane of OHS, which is the region where the OH 
airglow is suppressed. A total of 6 frames was acquired each day. FS 135 was observed after the 
observations on March 1 and 2, and FS 141 on May 18 for flux calibration. The seeing size was 
tf'7-tf'8. 

JH-band spectra were acquired using a C//95 slit, providing resolution of 210 at 1.65/um. Single 
exposure time was 1000s and the telescope was nodded after each exposure by 10" along the slit. 
The position angle of the slit was again set at —69° on March 1 and 2 to observe Jl and J2 
simultaneously. On May 18, the position angle was set at 24° to observe Jl and Jin. Total 
exposure time was 6000s on March 1, 4000s on March 2, and 8000s on May 18. F0 star SAO 
120154, AO star SAO 120250, and F5 star SAO 120186 were observed after the target on March 1, 
2, and May 18 respectively, for correcting atmospheric and instrumental transmission curve. 

The log of the spectroscopic observations are summarized in Table 1. 

2.3. Data Reductions 

The imaging data are reduced in a standard procedure of flat-fielding, sky subtraction, bad- 
pixel correction, residual sky subtraction, and shift-and-add. Photometry is performed in an aper- 
ture of same size as the slit width and length used in the JH-band spectroscopy. The position 
angle of the aperture is set at —69°. Measured flux is summarized in Table 2. 

The spectroscopic data are reduced in a standard procedure of flat-fielding, sky subtraction, 
bad-pixel correction, and residual sky subtraction. They are then corrected for the atmospheric 
and instrumental transmission curve, and are coadded to create final frames. Spatial widths of 
the aperture to extract 1-D spectra are 2 / /2 for Jl and l'/7 for J2. In the if -band, wavelength 
calibration is carried out using the airglow lines in the raw frames. In the JH-band, because 
almost all the airglow lines are eliminated significantly by OHS, wavelength calibration is done 
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using a pixel-wavelength relation obtained by previous observations of CISCO without OHS where 
the airglow lines are available. Systematic error of this pixel-wavelength relation is estimated to 
be less than 0.5 pixel (< 3A). Flux of the spectra is calibrated by the photometry of the H- and 
if '-band imaging in Table 2. 

2.4. HST WFPC2 Imaging 

We also obtain archived HST WFPC2 F702W data of Abell 1835 field (PI: Kneib, J.-P; PID: 
8249), to investigate morphology and optical magnitude. All the frames with 2500s exposure each 
are combined into a single frame after cosmic-ray rejection. Total exposure time is 7500s. 

3. Results 

The final JHK spectra of Jl and J2 extracted from the February and March data are shown 
in Figure 1. Although the wavelength resolution is lower (A/AA = 200 — 400) than that of Tecza et 
al. (2004) (A/AA ~ 1500), the quality of the continuum spectra is greater in the Jii-band, owe to 
the high sensitivity of OHS. In Jl, we detect strong Ha+[N n]AA6548, 6583 and relatively faint [O 
n]A3727 and H/3. No [O m]AA4959, 5007 is found. In 32, Ha, H/3, and [O II] are detected, but all 
the lines are weaker than those in Jl. No line is resolved with the low resolution of OHS/CISCO, 
which is 750 km s _1 in the iT-band and 1400 km s _1 in the .ff-band. The JH-b&nd spectrum of 
Jl taken in May is identical to that taken in March, with weak H/3 and [O n] without [O in] lines. 
Jin is not detected. 

We next carry out line fitting with Gaussian profiles to estimate their redshift, flux and equiv- 
alent width. We also deblend Ha and [N n], which are merged due to the low resolution, assuming 
[N li]A6583/[N n]A6548=3.0. Because the S/N of the spectra is low, line widths are fixed to the 
resolving power determined by the slit width. Results are given in Table 3, and fitted profiles of 
Jl are shown in Figure 2. We also estimate an la upper limit of the line flux of [O m]A5007 in Jl 
to be 2.0 x 10 -20 W m~ 2 . The lines of J2 are weaker, and the S/N is quite low. However, because 
the redshifts of all the lines match within la errors, we conclude these detection to be real. 

4. Discussion 

4.1. Emission Line Cloud of Jl 

A diagram of Ha/[N n]A6583 and [O iii]A5007/H/3 is shown in Figure 3. The line ratios indicate 
that the emission line cloud of Jl is ionized/excited by massive stars, not by a strong AGN, which 
is consistent with the X-ray observations showing no sign of a luminous AGN unless it is Compton 
thick (Fabian et al. 2000). 
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Balmer decrement is measured to be Ha/H/3=4.0±2.2. If we assume the intrinsic line ratio of a 
star-forming region to be 2.9 (Case B, T = 10 4 K; Osterbrock 1989) and adopt the SMC extinction 
curve (Pei 1992), dust extinction is estimated to be E(B — V) = 0.361q'3q. Due to the large errors 
in the line fluxes, the value is quite insecure. Nevertheless, this indicates that the H n region suffers 
extinction of Ay ~ 1. 

Lower limit of the metallicity can be assessed by 03iV2 or R23 estimator. The metallicity 
derived from 03N2 is 12 + log(0/H) > 8.7 (Pettini & Pagel 2004), while that from #23 is 
12 + log(0/H) > 8.9 (Kobulnicky, Kennicutt, & Pizagno 1999). These results are consistent with 
that of Tecza et al. (2004), suggesting that this system is well polluted by metals supplied from 
the starburst activities. 



4.2. Gravitational Lensing Model 

The HST-F702W image is shown in Figure 4. It can be seen that Jl consists of many sub- 
components. We name them to be Jla, Jib, Jlc, and Jld as labeled in the figure. 

Recently, Smail et al. (2003; private communication) found Ca H+K and other absorption 
lines at z = 0.25 in the optical spectrum of Jl, which suggests Jl to be contaminated and strongly 
lensed by a foreground galaxy. Also, as seen in the Ha image (Figure 1 of Tecza et al. (2004)), 
the line brightness becomes weak at the position of Jlc and a "hole" is seen there. These results 
indicate that Jlc is a member of Abell 1835 at z = 0.25 and causing a strong lensing effect on Jl 
at z = 2.565. 

Downes k, Solomon (2003) claimed that estimated brightness temperature of CO(3— 2) is much 
smaller than that expected from the CO line ratio, and argued that this is due to the gravitational 
lensing effect. They roughly estimate the magnification factor to be ~25. They also constructed a 
simple gravitational lensing model using the H ST-F702W image, assuming that Jla, Jib, Jin and 
J2 are all the lensed image of a single object. However, this model failed to reproduce the position 
of the observed components. In addition, it is known that the redshifts of Jl and J2 differ slightly 
(Ivison et al. 2000; Tecza et al. 2004) and they are thought to be different objects. We therefore 
propose a new gravitational lensing model assuming that Jla, Jib, and Jld are lensed images of a 
single object by Jlc and Abell 1835. Hereafter, we refer to the lensing galaxy at z = 0.25 as "Jlc" 
and distinguish it from the lensed galaxy Jl at z = 2.565. 

4.2.1. Lensing Model by GRA VLENS 

We use gravitational lensing software GRAVLENS (Keeton 2001). The model potential con- 
sists of two components Jlc and Abell 1835, both lie at z = 0.25. The potential of Jlc is assumed 
to be an isothermal ellipsoid with core radius 6 = 0, while that of Abell 1835 an NFW model 
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with parameters taken from Schmidt et al. (2001). We then carry out fitting calculation of the 
positions and aperture fluxes of Jla, Jib, and Jld, by varying velocity dispersion a, ellipticity e, 
and position angle of the ellipticity 9 e of Jlc. All the parameters of the lensing model is shown in 
Table 4. The uncertainty of each parameter is evaluated as follows: We start at the best fit model, 
and vary a single parameter at once. For the varied model, x 2 fitting is carried out by varying the 
position of the source to find the best-fit position. The average displacement between the positions 
of the best-fit spots and the original spots in the image plane is then calculated, and the error for 
the parameter is defined to be the value where the average displacement becmoes less than O'/l. 

Using this model, we also reconstruct the lensed image of Jin, J2, and the extended Ha cloud 
found by Tecza et al. (2004). Positions and magnification factors of these components are listed 
in Table 5. 

Figure 5 shows the image- and source-plane of the model. The solid contours show the profile 
of Jl at z=2.565 with an intrinsic diameter of O'/l (0.8 kpc). It can be seen that Jla and Jib in 
the model image are elongated, while that in the HST-F702W image is not. This suggests that 
the intrinsic size of Jl is smaller than O'/l (0.8 kpc), and that what we are seeing is a kind of a 
circumnuclear starburst disk like that in a nearby ULIRG. The total magnification factor of Jl is 
found to be 34^ . However, we should keep in mind that this magnification factor is that of a 
point source. If the source is assumed to be extended by O'/l, the total magnification factor becomes 
26±* 4 . 

Model contours of Jin, J2, and the Ha cloud are plotted with dotted lines. The extended 
Ha emission line reported in Tecza et al. (2004) is well reproduced by a cloud with an intrinsic 
size of 0'/2xO / /4 (1.6 kpc x 3.3 kpc) located at O'/l(0.8 kpc) west of Jl. Its total magnification factor 
estimated from our model is 23lg assuming that the Ha emission line region is extended by Of! 3. 
This value matches that of the CO emission which is estimated to be ~ 25 (Downes & Solomon 
2003). 

Jle, a counter image of Jl with half the flux of Jld, as well as Ha-d should appear in the 
west side of Jlc. However, both are not detected in either the HST-F702W or the Ha image. 
There are several explanations for this that; (1) they are absorbed by dust in the line of sight, 
presumably in Jlc, (2) the real potential of Jlc deviates from the model and the flux of Jle and 
Ha-d are depressed, and (3) Jla, Jib, and Jld are not a multiple image of a single source, but 
different components of a merging source structure. In the last case, the mass of Jlc should be 
much smaller to avoid the multiple lensing. However, because the optical properties of Jlc well 
matches our model as discussed in section 4.2.2, it is unlikely that the mass of the Jlc is much 
smaller than our model. We therefore assume that Jle and Ha-d are depressed either by dust or 
by shape of the real potential. 
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4-2.2. Properties of the Lensing Galaxy 



To check the feasibility of the gravitational lensing model, we next assess the optical properties 
of the lensing galaxy Jlc and compare them with the model. 

The profile of Jlc in the HST-F702W image is fitted with the Sersic (1968) profile 



where b n = 1.9992n — 0.3271 (Graham 2001). Fitted parameters are n = 2.1 and r e = 0?5, and 
the flux within r e is estimated to be 3.1 fiJj. Using spectra of 2-5 Gyr instantaneous burst model 
calculated by PEGASE (Fioc Sz Rocca-Volmerange 1997), the stellar mass is then estimated to be 
1.5-3.8xlO 9 M . 

On the other hand, the projected mass within r e of Jlc estimated from the gravitational lensing 
model is 3.8xlO 9 M using the equation of Kochanek (1995), which contains both stellar mass 
and dark matter halo. Considering the dark matter fraction in early-type galaxies of < 0.5 (Treu 
& Koopmans 2004) , this is a good agreement with the value obtained above and our gravitational 
lensing model seems reasonable. In addition, the size (r e =0'. / 5 ~ 2 kpc) and the velocity dispersion 
(a =123 km/s) indicate that Jlc is a relatively small elliptical in Abell 1835. 



To quantify the stellar contents of Jl and J2 together with their star formation history, we 
combine the optical BRI photometry of 100 with the current spectra to construct SEDs from the 
optical to NIR. Because 100 use 3"aperture for the photometry while our slit size is much smaller, 
we measured the ratio of the flux between the different aperture sizes in the .R-band, using the 
HST-F702W image convolved to the 0'/7 seeing size of the JH-b&nd spectroscopy. This ratio is 
used to correct all the fluxes of the optical bands in 100. Also, an additional error of 0.1 mag is 
added to the optical data considering the systematic uncertainties in the flux calibration of the 
NIR spectra. The emission lines in the NIR spectra are removed by the Gaussian fitting described 
in section 3. The obtained SEDs are shown in Figure 6. Continuum breaks at rest-frame 4000A 
are seen in both Jl and J2. 



We fit the SED of J2 by spectrum templates calculated by PEGASE (Fioc & Rocca-Volmerange 
1997) to which an empirical dust extinction law of starburst galaxies by Calzetti et al. (2000) is 
applied. Free parameters are dust extinction E(B — V), age, and total mass. The initial mass 
function of Salpeter is adopted. We consider two extreme cases of star-formation history; one is 




4.3. SEDs and Model Fitting 



4-3.1. Component J 2 
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an instantaneous burst, and the other is a continuous burst with duration of 1 Gyr. The results 
are summarized in Table 6, and the fitted spectra are shown in Figure 6b. J2 is expected to be 
dominated by very young population, and have formed of order of ~ W 8 M & stars. We must keep 
in mind that the SED is extracted from a small aperture (l"x2"), and what we see is only the 
central region of J2. 3"aperture flux in the i^-band (100) indicates that the obtained mass may be 
all underestimated by a factor of ~ 3. 

4-3.2. Component Jl 

To assess the contamination of Jlc to the SED of Jl, the spectral shape of Jlc is assumed to be 
that of a 2Gyr instantaneous burst model. It is scaled by the flux obtained from the HST-F702W 
image, and plotted in Figure 6a as the dotted line. This model spectrum of Jlc is then subtracted 
from the SED of Jl to obtain the intrinsic SED, which is shown as the dotted plots in Figure 6a. 

It can be seen that the U-band flux is little contaminated by Jlc. This agrees with the result 
of 100 that the peak of Jl in the [/-band is slightly offset toward east from that in the i?-band, 
indicating that Jl at z = 2.565 dominates the U- and 5-band flux. 

We next estimate the total intrinsic flux of Jl in the .fT-band to be 25 (iJy, by subtracting the 
model spectrum of Jlc obtained above from the K-band flux of 100 (3"aperture). Assuming the 
magnification factor of Jl to be 5 as a lower limit and the shape of the SED to be same as that 
of J2, the upper limit of the stellar mass in Jl is estimated to be ~ 10 9 M Q . If more dust exists in 
this system than in J2, this upper limit of the stellar mass will become larger. 

4-3.3. Formation Site of a Normal Galaxy? 

The amount of gas reservoir in Jl is fewxlO 9 M , estimated from the CO observations (Downes 
& Solomon 2003). Therefore, if SMM J14011+0252 is a place where a massive galaxy of ~ 10 12 M Q 
is being assembled, stellar mass of ~ 10 12 Mq must have already been assembled and hidden in 
dust. This is more than 100 times as much as the value calculated above, and extinction of 
Ay > 5 is required among the galaxy. Although such an amount of extinction is broadly observed 
in the circumnuclear starburst disk in nearby ULIRGs (e.g. Genzel et al. 1998), their size is 
smaller (< 0.5kpc) and their stellar and gas masses are comparable to each other and much smaller 
(~ 10 10 Mq) than those required in SMM J14011+0252. Therefore, considering the difficulty to 
hide a stellar component of 10 12 M & in dust, it is more probable that SMM J14011+0252 is a site 
of less massive (< 10 10 Mq) galaxy formation. 

If we assume that the submillimeter radiation is emitted from the same region that emits CO 
or Ha and its magnification factor is of order of 30, the intrinsic 850 /jm flux is ~ 0.5mJy. This 
is the flux that most of less massive high-z galaxies such as Lyman break galaxies are expected 



-10- 



to emit at most (Chapman et al. 2000) . Also, the luminosity in the far-IR wavelength becomes 
~ 10 12 Lq, which are same as that of local ULIRGs. 

However, readers should keep in mind that above discussions are based on the uncertain 
gravitational lensing model, which does not have redshift information of the observed components 
to settle the membership. The uncertainty in the model parameters are also large and the calculated 
total magnification may change by factor of 2. Also, there is a possibility that that the position of 
the submillimeter source is offset from that of CO or Ha and it does not suffer a strong lensing 
effect. In such case, the total magnification is only a factor of 5, and the intrinsic luminosity in 
far-IR becomes an order of magnitude larger. Then, this object may be interpreted as a massive 
and dusty starburst forming a larger mass (~ 10 11 M Q ) of stars. 

It is likely that the starburst in this system are triggered by an encounter of Jl and 32. If we 
ignore the radial distance (that is, assuming the distance between Jl and J2 is that observed in the 
source plane of l'.'2=10 kpc) and assume the encounter to have occurred 50 Myr ago, the projected 
relative velocity is estimated to be 200 km s _1 . Considering the radial distance and motion of Jl 
and J2, this matches the typical relative velocity of 1000 km s _1 in a cluster of galaxies. 

Therefore, Jl may be a high-z counterpart of the local ULIRGs such as Arp 220, and its 
submillimeter flux are magnified above the detection limit of SCUBA by the strong lensing effect 
of Jlc. Its total barionic mass is ~ 10 10 M , and we hypothesize that this object is a less-massive 
forming galaxy at z ~ 3 which will evolve into a mid-sized elliptical or spiral galaxy such as the 
Milky-Way. 

5. Summary 

We have conducted JHK spectroscopy of components Jl and J2 of the gravitationally lensed 
submillimeter galaxy SMM J14011+0252 at z = 2.565 using OHS/CISCO at the Subaru telescope. 

Ha, H/3, [O n] lines in both Jl and J2 are detected , and [N n]AA6548, 6583 lines are also 
detected in Jl. No [O in] lines are found. From the emission line diagnosis, the emission line cloud 
of Jl is shown to be ionized/excited by stellar UV radiation, not by an AGN. The continuum is 
detected in both Jl and J2, and the Balmer breaks are clearly seen. The SED of J2 is fitted by a 
model spectrum of young (~ 50 Myr) stellar population with some amount of dust extinction. 

The new gravitational lensing model is constructed to reproduce the small knots of Jl, which 
also reproduces the morphology of the extended Ha cloud well. This model suggests that Jl may 
be strongly magnified by the elliptical galaxy in Abell 1835 by a factor of 30. Also, using this 
model together with the profile fitting of the lensing galaxy, the stellar mass of Jl at z = 2.565 
is estimated to be < 10 9 Mq. This indicates that SMM J14011+0252 is not a formation site of 
a present-day massive elliptical, but more like a high-z version of an ULIRG, and will become a 
present day spiral or mid-sized elliptical of 1O 1O_11 M0. 
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Fig. 1— JHK spectra of (a) SMM J14011+0252 Jl and (b) J2. Both spectra are smoothed by 3 
pixels bin. Dotted lines indicate 1 a error level. 
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Fig. 2.- 
lines) . 



(a) Ha+[N n], (b) H/3, and (c) [O n] lines of Jl, fitted with Gaussian profiles (dashed 
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Fig. 3. — Emission-line diagnostic diagram of Ha/[N n]A6583 vs. [O iii]A5007/H/3. The large open 
circle with error bars is the line ratio of Jl, while the large filled circle indicates that of SMM 
J04431+0210 taken from Prayer et al. (2003). The dotted lines show the separation between H 
II regions, LINERs, and Type 2 AGNs taken from Vcilleux & Osterbrock (1987) and Filippenko 
& Terlevich (1992). Also the line ratios of Seyfert 2s(filled circles), radio galaxies (filled squares), 
narrow emission-line galaxies (asterisks), H n regions(open circles), H n galaxies(open circles), and 
starburst galaxies(open triangles) taken from literatures (Veilleux & Osterbrock 1987; McCall et 
al. 1985; French 1980; Koski 1978; Costero & Osterbrock 1977) are plotted with smaller marks. 
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Fig. 4.— HST-F702W image of SMM J14011+0252. The dotted line indicates the slit position of 
the OHS/CISCO spectroscopy. 
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Fig. 5. — The gravitational lensing model image of (a) the image plane overlaid on the HST- 
F702W image and (b) the source plane. The solid contours show the profile of Jl at z=2.565 with 
an intrinsic diameter of O'/l (0.8 kpc), while the dotted contours show that of Jin, J2, and the 
Ha emission line cloud. Intrinsic sizes of Jin and J2 are Of! 2 (1.6 kpc), while that of the Ha cloud 
is an ellipse of / .'2xO / . / 4 (1.6 kpc x 3.3 kpc). 
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Fig. 6. — Optical to NIR SEDs of (a) Jl and (b) J2. Present data are shown as open circles, while 
the data taken from 100 are shown as filled circles, (a) The dotted line shows the assumed model 
spectrum of z = 0.25 lensing galaxy, and the dotted plots show the lens-subtracted SED. (b) The 
dashed and dotted lines show the result of the fitting by instantaneous burst and 1 Gyr continuous 
burst models, respectively. 
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Target 


Date 


Band 


texp 


PA 








(s) 


(°) 


J1/J2 


2002 Feb 28 


K 


4800 


-69 


J1/J2 


2002 Mar 1 


JH 


6000 


-69 


J1/J2 


2002 Mar 2 


JH 


4000 


-69 


Jl/Jln 


2002 May 18 


JH 


8000 


24 



Table 1: Log of spectroscopic observations. 



Component 


Aperture 


F x ,k> 


F X ,H 






(xlO" 18 W m" 2 /an" 1 ) 


(xl0~ 18 W m" 2 fim.- 1 ) 


Jl 


0"95x2"2 


9.2±0.5 


12.2±0.7 


J2 


0795 xl?7 


2.4±0.4 


3.8±0.5 



Table 2: Flux in the K' and #-band of Jl and J2. 





Line 


^rest 


Line Flux 


Continuum Flux 


EW rcst 


z 






(A) 


(xl0~ 19 W m~ 2 ) 


(xlO" 18 W m~ 2 ^m" 1 ) 


(A) 




Jl 


Ha 


6563 


1.3±0.4 


8.9±1.8 


42 


2.565±0.001 




[N II] 


6583 


0.62±0.39 


8.9±1.8 


20 


2.565±0.001 




H/3 


4861 


0.33±0.12 


10.9±0.5 


8.5 


2.562±0.002 




[0 ii] 


3727 


0.41±0.13 


14.0±0.6 


8.1 


2.563±0.003 


J2 


Ha 


6563 


0.48±0.34 


1.9±1.1 


72 


2.562±0.003 




[N ii] 


6583 


0.27±0.33 


1.9±1.1 


40 


2.562±0.003 




H/3 


4861 


0.11±0.09 


3.5±0.4 


9.0 


2.561±0.005 




[0 ii] 


3727 


0.26±0.09 


3.6±0.4 


20 


2.564±0.003 



Table 3: Fitted parameters of the emission lines in Jl and J2. 



Component Parameter Comment 

Jlc a 123^4 km/s Velocity dispersion of the potential 

e 0.671q;5o Ellipticity of the potential 

#e ^-io ° Position angle of ellipticity 

b 0+ - 07 " Core radius 

Abell 1835 7 S 0.64+f Mpc 

k s u.ioy_ 0015 

dx 44^4 " Distance from Jlc in R.A. 

dy 20^4 " Distance from Jlc in Dec. 



Table 4: Model parameters of the lensing potential. 
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Source Name 


Position 3, 


Image Name 


Magnification 


Position b 


Jl 


(24.71, 10.64) 


Jla 


11 


(-0.59, -0.18) 






Jib 


13 


(-0.25, - 0.67) 






Jld 


6.3 


(-0.23, 1.10) 






(Jle) 


3.1 


(0.46, -0.21) 


Jin 


(24.38, 11.10) 


Jin 


4.4 


(-1.00, 2.25) 


J2 


(25.79, 10.93) 


J2 


5.7 


(1.89, 0.80) 


Ha 


(24.62, 10.57) 


Ha-a 


12 


(-0.78, 0.28) 






Ha-b 


10 


(-0.44, -0.73) 






Ha-c 


9.5 


(-0.51, 0.90) 






(Ha-d) 


2.0 


(0.36, -0.06) 



Table 5: Calculated parameters of the lensed images. The position of each components is shown 
in Figure 5. The parenthesized images are those not detected. Note that these magnification 
factors are calculated asuuming that a lensed image is a point source. If the image is extended, the 
magnification factor becomes smaller (see text). 



"Position on the source plane in arcsec, relative to Jlc. 
^Position on the image plane in arcsec, relative to Jlc. 



Model 


Age 


E(B - V) 


Stellar Mass 


Total Mass 


SFR 


x 2 




(Myr) 




(M e ) 


(M e ) 


(M Q yr- 1 ) 




Instantaneous 


8 


0.33 


6.1xl0 7 


6.2xl0 7 


0.0 


16.85 


lGyr Continuous 


70 


0.40 


1.7xl0 8 


2.4xl0 9 


4.4 


22.83 



Table 6: Parameters for the best-fit models of the J2 SED. Mass and SFRs are corrected for the 
gravitational lensing effect with the magnification factor of 5.7. 



